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Abstract: 

We present an exact a s calculation of the Wilson coefficients associated with the 
dipole moment operators. We also give an estimate of the branching ratio for b — > 57. 
We find that higher dimensional effects are under control within 9% for BR{b — > 57) = 
(4.3 ±0.37) x 10~ 4 . 
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In an article 1 by us recently, we presented a leading logarithmic analysis of the heavy 
particle effects on the process b — > S7, which incorporates a complete operator mixing. 
Among the results, we find that our mixing matrix differs in some elements from the 
work by Misiak 2 . Had we subscribed to the prescription of further mixing with evanescent 
operators 3 , we would have obtained extra contributions (in unit of fir 2 /87r 4 ) 

A 70 67 Obi = 1 A A 70 67 O s2 = -8/3, 

A 70 68 Obi = -1/2, A 7 o 68 52 = -4. (1) 

These would have made the mixing by Ciuchini et. al. 4 coincide with ours. (We thank 
Misiak for correspondence on this.) For earlier work, see 5 ' 6 . 

More importantly, by choosing two different limits for extrapolation: rat = rn w and 
m t 3> m w , we estimated the effects due to higher orders in m 2 /m 2 to be about 20%. In 
view of recent interests in the experimental branching ratio 7 for B s — > K* + 7 and an 
impending value for the inclusive rate B s — > X s + 7, together with the sensitivity of these 
processes as a short distance probe, an uncertainty in short distance analysis of 20% is 
hardly satisfactory or even acceptable. The purpose of this publication is to provide some 
remedy. 

One can trace the uncertainty to the fact that m 2 w jm 2 t « 25% for m t ~ 150 Gev. It 
strongly suggests that one should calculate the order a s diagrams exactly, whereever there 
are internal top and/or W-boson propagators. This is what we have done. 

In all two-loop diagrams (and their attendant counter terms) which contribute to 
this calculation, we keep all orders in x = m^/m 2 ^ but discard terms proportional to 
0{m\/m1 w ) or 0(m 2 /m 2 w ) (we have factored out the Fermi weak coupling constant). 
The amount of algebra is highly non-trivial and is aided by Schoonschip. 8 
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Given a diagram in which there are top and/or W — internal lines, there are various 
sequences of operations one can follow to isolate its dependence on the heavy masses 9 . In 
any case, the underlying method is based on partitioning of the diagrams into heavy parts 
and operator inserted matrix elements. For the present situation, we treat the top and the 
W — 4> as being correspondingly heavy, relative to other masses and external momenta. 
In the order we are working with, a heavy part always contains both the top and the 
W — <fi internal lines. Vertices made of light particles and momenta acting on them and 
Wilson coefficients which contain all dependence of heavy masses are organized within this 
formalism. In this way, we obtain the effective Lagrangian 

Leff = J2 C ^ ( 2 ) 

i 

where Cj are the Wilson coefficients and Oi are a set of local operators, made of light 
fields. 

Of particular interest in the present order a w a s calculation are the coefficients Co 51 
and Co 52 , with the accompanying operators 

051 = ig s sG^a^ u {m s P L + m b P R )b/2, (3) 

and 52 = 51 (g s G^ -> (-l/3)eF^). 

Before we go on further, let us define 

Co 5l = W^Co^/Gt (i = 1,2), and G t = 2V2G F V*V tb , (4) 
where Gf is the Fermi weak coupling constant and Vs are the Cabbibo-Kobayashi- 
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Maskawa matrix elements. Then the exact a s results are (in the R* scheme to be defined) 
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where the last terms in big parenthesis in each of the two equations above are the contri- 
butions due to u and c quarks. Sp stands for Spence function. Eq.(5) are to be contrasted 
with the asymptotic results (m t ^> m w ) 
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and 

4 >, («ymW4a = -- + - + - log \-± j + - log ^ j , (6) 

obtained by discarding all 0(m^/m|). One can use Eq.(5) to check some of the mixing 
matrix elements of the renormalization group equations (RGE). Unfortunately, 7o 67 o 6 i> 
7o 67 o 52 5 7o 68 o 5 n and 7o 68 o 52 do not enter to this order. Their determination has to 
come from three loop diagrams. In our opinion, direct Feynman diagram computation of 
Green's functions for processes to extract mixing matrix elements should be the definitive 
procedure. 

We plot the exact a s result Eq.(5) in Figure 1, together with the asymptotic result 
Eq.(6) for Cq^ 2 % /4a s . For x=4, the discrepancy is in fact about 50%. It is interesting to 
note that the exact Cq^ is quite flat between x=l to x=10. The dotted line is the a° s 
exact result, also known as the Inami-Lim functions 10 

=( ) _ -2x - 3x 2 + 6x 3 - x 4 - 6x 2 logx 
°" ~ 4(1 -x) 4 ' 



- (o) _ 7x - 12x 2 - 3x 3 + 8x 4 + (12x 2 - 18x 3 ) log x 
C °*» ~ 4(1 - xY ' ( } 



We see that below x=6.9, second order QCD correction is bigger than the lowest order 
result. This has been known for the approximate results for some time and in fact is an 
impetus for looking into rare decays of this genre. 

Please note that if one uses Eq.(7) as (a part of) the boundary conditions, then Eq.(5) 
are the a s solution of RGE for all values of x, in so far as the overall coefficients to log/U 2 



are concerned. We would like to stress that we are using the R* scheme. In the MS 
scheme, the exact result to order a s is obtained if one makes the replacement m t — > 
m t (l — a s /ir log(m 2 /// 2 )) in Eq. (7); this replacement is just a finite renormalization to 



go from the R* scheme to the MS scheme. 
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We shall make the assumption that after a s corrections, it is safe to add the leading 
logarithmic terms to complete the leading QCD sum. In other words, we assume that the 
higher order QCD corrections can be obtained either in the limit m^/m 2 <C 1 or in the 
limit m 2 /m 2 = 1. This assumption can be tested as in our previous publication. For 
(7o S2 , there entail two different extrapolations 



(exact) + C n Q l ner ° raer (m t = m w ) (8.a), 

and 

4°) 2 +(1) (exact) + h o l her order (m t » m w \ (8.6) 

where er order (m t = m w ) and er order (m t > m w ) are the remaining leading log- 
arithmic sums with the boundary conditions set at m t = m w and m t ^> m w , respectively. 1 

For the physical process b — > s+7, some four quark operators also contribute, resulting 
in an effective coupling 2 

1 3 

= C 5 2 + -^2 C 067 + -^2 C 068 ■ ( 9 ) 

Also, to remove the dependence on |Gt| 2 , which is not accurately known experimentally, 
we normalize the b — > S7 partial width to the well established semileptonic b — > cev partial 
width, and use the following relation 11 |V^Vth| — |Vc&|. 
This ratio is given as: 

Lj^Jt K a -9SR ( 1 - ^pl f(m j mh )) " ic&W, (io) 

r (b — > cev) 6it g (m c /mb) \ on J ° 2 

where g(m c /mb) ~ 0.45 and f(m c /mb) — 2.4 correspond to the phase space factor and 
the one-loop QCD corrections to the semileptonic decay, respectively. 

In Figure 2, we have plotted BR(b — > S7) as a function of mt. The solid and dashed 
curves are obtained with the aid of the interpolation equations in Eq.(8), together with: 

BR {b _» S7) = L£ n s j1 br {b _> ceV) (ii) 

1 (0 — > cev) 
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BR(b^ceu) ~ 0.108. (12) 

The vertical dotted line to guide the eyes intersects these curves at m t = 140 GeV and 
gives respectively: 



BR(b — > S7) =4.66 x 10" 4 , 

BR(b — > S7) = 3.93 x 10" 4 . 
From Eq. (13), we obtain the following mean value (this is our estimate), 



(13) 



Si? (6 -> 57) = (4.3 ± 0.37) x 10" 4 . (14) 

This is to be compared with an upper limit 5.4 x 10~ 4 given recently by the CLEO 
Collaboration. ( 7 ) 

The uncertainty due to subleading logarithmic and higher dimensional effects is about 
9%, which is a big improvement and more reliable over what we gave before, where m 2 /m 2 u 
effects at or} ' were not treated. 

We now give some technical details. We use the general linear covariant gauges 
^(d^G^) 2 for the gluons. The complete cancellation of a for Co 51 and Co 52 is a stringent 
confirmation on the correctness of the algebra. The gauge fixing for W-fields is —C + C~ 
with C + = —d^W^ + m w (f) + + ieA^W^. For oversubtractons and renormalization, we 
use the i?*-scheme. Thus, let V be a one light particle irreducible (1LPI) diagram which 
contains the heavy particles, and let 7 represent a 1LPI graph or subgraph (7 C T) with 
external generic momentum p. We define 

T 7 = P°l e part 0/7, e = n — 4, 

I o 777, 

} = 1(P = 0) + p- 7 (p = 0) + • • • + - y P m Q^l{p = 0). (11) 
The R* renormalization procedure is defined as 

R*{lheavy) = (1 ~ T^ m) )-f heavy , R* ("flight) = (1 ~ flight, (12) 



where m is so chosen that the neglected terms are genuinely of 0(l/m\ e ). It is important 
to repeat that the results Eqs.(5,6) are given in the renormalized top mass under the R*- 
scheme, where j3 mt = 0. 



Except for trivial factorizable cases, all two loop integrals we need are related to 12 
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h,i,i(ml,ml,ml) = J d n kd n q 
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= ^(1 " a + b) 2 - 4b. 



(13) 



(14) 



(15) 



Assuming a, b ^> 1, (i. e. we take rri2 = m t , ms = m w ^> mi = m^), we expand 
/(a, 6) in series of m 2 , which must be retained to proper orders in intermediate steps. 

Details of this work are to be published 1 . 
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FIGURE CAPTION: 



Figure 1: Co 52 ( m b) dependence on nit to order a s with mj, = 4.8 GeV ', m TO = 81 GeV, 
oe s (m,b) = 0.19. See text for explanation of various curves. 

Figure 2: Branching ratio for b — > S7 as a function of m 4 with QCD corrections. The 
solid line represents the interpolation given by Eq. (8. a); the dashed line represents the 
interpolation given by Eq. (8.6); the dotted line represents the values of the branching 
ratio for m t = 140 GeV. We used m b = 4.8 GeV, m w = 81 GeV, a s (m b ) = 0.19. 
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